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Abstract 



We construct a new many-body Hamiltonian with two- and three-body 
interactions in two space dimensions and obtain its exact many-body ground 
state for an arbitrary number of particles. This ground state has a novel pair- 
wise correlation. A class of exact solutions for the excited states is also found. 
These excited states display an energy spectrum similar to the Calogero- 
Sutherland model in one dimension. The model reduces to an analog of the 
well-known trigonometric Sutherland model when projected on to a circular 
ring. 



PACS: 03.65.Ge, 05.30.-d 



Typeset using REVTp^X 



*Permanent Address: The Institute of Mathematical Sciences, Madras 600 113, India 
^Permanent Address: Indian Institute of Science, Bangalore 560012, India 

1 



There has been a renewed interest in exactly solvable many-body problems in recent 
times. A celeberated example of a solvable many-body problem is the well-known Calogero- 
Sutherland model (CSM) in one dimension |]]-f§. The particles in the CSM are confined in a 
one-body oscillator potential or on the rim of a circle, and interact with each other through a 
two-body potential which varies as the inverse-square of the distance between particles. The 
CSM and its variants in one-dimension, like the Haldane-Shastry model for spin chains [f§], 
have provided the paradigms to analyze more complicated interacting systems. It is therefore 
of great interest to find analogous models in higher dimensions which share features similar 
to the CSM. A characteristic feature of the CSM is the structure of the highly correlated 
wave function. The correlations are built into the exact wave function through a Jastrow 
factor ( \ a for any pair of particles denoted by The exponents on the 

correlator are related to the strength of the inverse-square interaction. Notice that this 
factor is antisymmetric (symmetric) in particle labels for A = 1(0) and vanishes as the two 
particles approach each other. A generalization of this in two dimensions is to be found in 
the Laughlin's trial wave function where the correlations are built in through the factor 
(zi — Zj), where Zi are the particle coordinates in complex notation. The corresponding 
Hamiltonian for which the Laughlin wave function is an exact ground state has not been 
analyzed to the same degree of detail as the CSM. However it is known that it is the ground 
state for a Hamiltonian describing spin polarized electrons in the lowest Landau level with 
a short-range repulsive interaction [[J. 

In two dimensions, however, there exists another form of the pair correlator with which 
a Jastrow-type many-body wave function may be constructed, namely 

Xij = Xii/j - xjUi . (1) 

This type of correlation is by definition antisymmetric and goes to zero as two particles 
approach each other. In addition, it introduces zeros in the wave function whenever the 
relative angle between the two particles goes to zero or tt. The difference with the Jastrow- 
Laughlin form is also significant; ([]]) is a pseudoscalar. Unlike the Laughlin type correlation, 
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it does not impart any angular momentum to the two dimensional wavefunction. We now 
ask: Is there a Hamiltonian for which the exact solutions have this type of correlation? In 
this paper, we construct a many-body Hamiltonian whose exact ground state involves a 
product of the pairwise correlations of the form given by Eq. ([I]). While we are able to find 
a class of excited states of the Hamiltonian, we have not been able to solve for all the excited 
states. The known spectrum of states bears a remarkable similarity to the spectrum of the 
CSM in one dimension; namely, the spectrum (which includes the ground state and part 
of the excited states) is bodily shifted by an iV-dependent but configuration independent 
factor. This is the feature which is responsible for ideal fractional exclusion statistics |7||| 
in the CSM [|. 

The N— particle Hamiltonian which we consider is given by 

fi 2 " „ ~ . ft 2 " ff ft 2 ' fi.n 

1 = 1 1=1 l,J %J !J,t l J lK 

where is given by (jl]); g\ and g 2 are dimensionless coupling strengths of the two- and 
three-body interactions respectively. While g\ and g 2 can be independent of each other in 
general, for the type of solutions involving the correlator in (|l|) they are not. We will specify 
their relationship shortly. The particles are confined in a one-body oscillator confinement 
potential. The Hamiltonian is manifestly rotationally invariant and symmetric in all particle 
indices. As in the CSM, we may scale away the mass m and oscillator frequency uj by scaling 
all distances r% — > ^Jmu/h fi, and measuring the energy in units of Too. This is done by 
simply setting H = m = u = 1. In these units the Hamiltonian is given by 

1 N 1 N N f 1 N ff ff 

1 i=l 1 i=l 1 i,j X ij 1 i,j,k ^ij^ik 

We first obtain the exact many-body ground state of this Hamiltonian. As an ansatz for 
the ground state wave function, consider a solution of the form 

N i iV 

M*i,Vi) = II X ?j \ X a\ a exp(--^rf) . (4) 

i<j i=l 

The ansatz for the ground state is bosonic (fermionic) for A = 0(1). Let us define 
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g = a + A . (5) 

Clearly \l/o correctly incorporates the behavior of the true ground state in the asymptotic 

region, \fi\ — > oo and for g > 0, \l/o is regular. The eigenvalue equation now takes the form 

1 N r 2 -- 1 N r-fu 

H* = [-( gi -g(g-l)) ]T ^ + -{g2~g 2 ) £ + ^(iV - 1) + iV]* . (6) 

Therefore \l/o is the exact many-body ground state for arbitrary number of particles of the 
Hamiltonian if 

gi = g(g-i), 92 = g 2 - (7) 

Equivalent ly, 



which immediately implies g\ > —1/4 for physical states and #2 is positive definite. Note 
that these relations are identical to the ones obtained in the CSM. While \l/o is a solution 
of the Schrodinger equation in either of the two branches given in (||), the solutions in the 
lower branch are regular only in the limited range — 1/4 < g\ < 0. The solutions in the 
upper branch,however, are regular for all g\ > —1/4. With this the ground state energy is 
given by 

E = N + gN(N-l). (9) 

Note that this has exactly the form of the ground state energy of the CSM. 

We emphasise that our objective here is not to find the general solutions for arbitrary 
gi and g2, but to find a Hamiltonian whose solutions have the novel correlation in Eq. ([!]) 
built in. With the form of g\ and #2 given above in (f?]), the solutions found above are indeed 



the lowest energy states. To see this [ 10 1 , let 



y ■ 

Q Xl = p Xl - ixi + i g 22 Y~ » 

iU/i) ij 

Qvi = Pvi - Wi - ig E -y~ ■ ( 10 ) 



It is easy to see that both Q Xi and Q yj annihilate the ground state. The Hamiltonian can 
be recast in terms of these operators as 

\ E [QlQ*i + QlQvi = h - Eo, (ii) 

where E is given by Eq. (||). Clearly the operator on the left hand side is positive definite 
and annihilates the ground state wave function given by Eq. (|4]). Therefore E must be 
the minimum energy that an eigenstate can have. Note that the total angular momentum 
L = J2A x iPyi — ViPxi) commutes with the Hamiltonian, and may be used to label the states. 

While we have not been able to find the complete excited state spectrum of the model, 
the eigenvalue equation for a general excited state may be obtained as follows. From the 
asymptotic properties of the solutions of the Hamiltonian in Eq. (|3p, it is clear that \I> has 
the general structure 

tf(a*,2/i) = V (xi,yi) $(xi,yi) , (12) 

where \l/o is the ground state wave function, and $ is a completely symmetric function of 
the coordinates. (Note that the bosonic or fermionic nature of the eigenstates is completely 
specified by the choice of A). $ satisfies the eigenvalue equation 

1 N N 1 Q Q 

[ "2 EV, 2 + E^-V, + 9 )]* = {E-Eo)*. (13) 

i=l i=l i,j l 3 & % 1 

(¥i) 

We now consider a class of excited states built on this ground state, namely the radial 
excitations. These states are of the form 

N 

*{{x uyi }) = y ({xi,Vi}Mty, * = E r * ( 14 ) 

i=i 

The eigenvalue equation for \1/ can then be easily reduced to a confluent hypergeometric 
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( _ii + [6 _ i] _Li_ aO(i) = 0. (15) 

where, b = E = N + (A + a)N(N — 1), a = —\[E — E ]. The admissible solutions 
are the regular confluent hypergeometric functions, = M(a,b,t). The normalizability 



immediately imposes the restriction a = —n, where n is a positive integer and $(t) is 
in general a polynomial of degree n. The corresponding eigenvalues follow immediately: 
E = Eq + In. Once again these solutions have the spectrum that is identical to that of 
states in the CSM. The tower structure of eigenvalues, separated by a spacing of 2, is also 
reminiscent of what happens in the case of anyons |JL2{| . Notice that until now we have only 
considered the completely symmetric solutions modulo the correlation factor Xy. Therefore 
the bosonic or fermionic nature of the eigenstates is completely specified by the choice of A. 

We now seek exact solutions to Eq. (fl3|) which are not in the same class as the simple 
solutions obtained above. We do this by eliminating the cross-terms involving the derivatives 
in Eq. ([13]). As in the case of the CSM ||, we assume that such eigenstates have the form, 

${Xi,yi) = a m 1 ...m N ,n 1 ...n N ]J H mi (Xi)H rH {yi) = ^ a ™i,ni II H ™A X i) H nAVi) ) 

rrn,n% rrii,ni rrii,ni rrii,ni 

(16) 

where m^rii are positive integers with the restriction J2A m i + n i) — N, and the H m 's are 
Hermite polynomials of degree m. These are of course solutions of the non-interacting N- 
oscillator Hamiltonian. The expansion coefficients a mitni are determined by demanding that 
the cross-terms vanish when acting on the wavefunction above ]TB[. Indeed the cross-terms 
vanish, provided 

4 {iTLi ~\~ l)(^j ~\~ l)Qmi...mi+i...m 3 '...rij...rij-|-i... 

"I" (jTT'i ~\~ l)<2 mi ,..m i+1 ...raj ...rn...Tij—\ ... (j^j l)^mi...mj— \...mj...rii...rij-\-\... 
4 [TTlj ~\- l)(?T.j -|- l)(3 mi ...m i ...m.jj r \ ...rii+\ ...rtj ... 

-\- {jUj ~~\~ .. n%i.. .mj+\...ni—\...nj... ij^i ~^)Q"m\...m,i...m,j^\...ni+i...nj ... i (^^) 

for all < rrii, < M . In general, it is difficult to obtain a closed-form solution for the 
coefficients a muni by solving the above equation. However, in specific cases involving the low 
energy excitations, we have verified that it is possible to satisfy the above identity and the 
solutions exist. It is not apparent, however, that such solutions exist in general. If they do, 
then the energy spectrum of these excited states is given by = i?o+X]i( m i+ n i) = E +J\f, 



where M now labels the quanta of oscillator excitations above the ground state. 

Finally, we point out the connection between our model Hamiltonian and the trigono- 
metric Sutherland model 0. Restricting the particles to move along the perimeter of a unit 

circle in the Hamiltonian (|3]) without the confinement potential, we get 
1 N pa N i N 

H = + + f E [!+«*<«, -»,)«*<«, -ft)], (18) 

(¥i) (¥i#0 
since Xij = sin(#j — 9j) now. Using the identity 

f] cot(g i - gj) cot(g t - fij] = - N ^ N ~ 3 )(iV ~ 2) , (19) 

i,j,k 

we immediately recover an analog of the trigonometric Sutherland model, but shifted by the 
constant g 2 N(N — 1)(N — 2)/3. Note, however, that the potential in flI8|) depends on the 
function sin(#j — 0^), rather than the chord-length which is proportional to sin[(0j — 9j)/2] . 

In summary, we have constructed a nontrivial many-body Hamiltonian in two dimensions 
with two- and three-body interactions which can be solved exactly for the ground state and a 
class of excited states. The wave function involves a novel type of correlation which vanishes 
not only when the particle positions coincide, but also when the relative angle is a multiple of 
7T. While the later may seem unphysical, the paucity of exactly solvable models even for the 
ground state in higher dimensions makes this an interesting system to study. The spectrum 
displays a dramatic similarity to that of the Calogero-Sutherland model in one dimension, 
namely, the ground state and the known excited states are all shifted by exactly the same 
quantity proportional to iV(iV — 1), where N is the number of particles. We also note that 
there are many similarities between the spectrum of our Hamiltonian and the spectrum of 
the many anyon Hamiltonian without, however, the complication introduced by the 



time-reversal violation in anyon physics. 
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